Incommensurate short-range order in the S=l triangular antiferromagnet NiGa2S4 
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Neutron scattering is used to investigate spin correlations in ultra pure single crystals of the S=l 
triangular lattice NiGa2S4. Despite a Curie- Weiss temperature of 0cw=-8O(2) K, static (t > 1 ns) 
short range (^a6=26(3) A) incommensurate order prevails for T>1.5 K. The incommensurate modu- 
lation Qo = (0. 155(3), 0. 155(3), 0), Qcw, and the spin wave velocity (c=4400 m/s) can be accounted 
for by antiferromagnetic third-nearest neighbor interactions J3=2.8(6) meV and ferromagnetic near- 
est neighbor coupling Ji =-0.35(9) J3. Inter-plane correlations are limited to nearest neighbors and 
weakened by an in-plane field. These observations show that the short range ordered glassy phase 
that has been observed in a number of highly degenerate systems, can persist near the clean limit. 

PACS numbers: 74.72.-h, 75.25.+Z, 75.40.Gb 
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In geometrically frustrated magnets no spin state opti- 
mizes all interactions and consequently spins can remain 
disordered at temperatures (T) well below the Curie- 
Weiss temperature, [ij The frustration results in high 
degeneracy as in strongly correlated electronic systems 
but in a context that is more amenable to analysis and 
comparison between theory and experiments. Frustrated 
magnets thus provides a unique opportunity to explore 
the new states of matter that can result from high de- 
generacy and suppression of conventional mean field be- 
haviors. The simplest frustrated magnet consists of an- 
tiferromagnetically interacting spins on a triangular lat- 
tice. For S=l/2 quantum fluctuations play an impor- 
tant role with long-range order only predicted to exist 
at T=0. [2-4] In addition for weakly distorted S=l/2 
systems close to a metal insulator transition such as n- 
(BEDT-TTF)2Cu2(CN)3, there is evidence of a spin hq- 
uid phase with no sublattice magnetization. 5] For the 
classical {S — 00) limit long range spin order at fi- 
nite temperatures has been observed in materials such 
as RbFe(Mo04)2 and NaMn02. [1,0] The situation for 
S=l is less clear with a resonating valence bond state and 
a spin nematic phase having been proposed. [8-10] Inter- 
estingly, there are no experimental examples of isotropic 
triangular lattice antiferromagnets with S < 1 which pos- 
sess conventional long-range order. It is therefore impor- 
tant to investigate the magnetic properties in isotropic, 
stoichiometric S=l triangular lattice systems, where col- 
lective phases may be possible. 

NiGa2S4 consists of two-dimensional (2D) planes of 
S—1 Ni^+ magnetic ions occupying the corners of 
edge sharing equilaterial triangles with Van der Waals 
forces alone stabilizing the lamellar structure. [Ill [l2| 
While magnetization measurements yield 6cw=-80(2) 
K, short-range (SR) order prevails for T>1.5 K. Here we 
study the low-T magnetic properties of NiGa2S4 using 



neutron scattering from high quality single crystals and 
provide clear evidence for quasi 2D incommensurate SR 
order. While proposals for a spin nematic phase may yet 
be relevant, we show that important aspects of the low 
energy physics are well described by a Ji- J3 model with 
strong third nearest neighbor (NN) antiferromagnetic ex- 
change and a weakly ferromagnetic NN exchange. In this 
model, the lack of long range would be associated with 
anomalous quasi-2D domain dynamics. 

The experiments used the SPINS and BT2 triple-axis 
spectrometers and the High Flux Backscattering (HFBS) 
instrument at the NIST Center for Neutron Research. 
Sintered pellets and two different single crystal config- 
urations were used. We used 7 co-aligned single crys- 
tals (with a mass of 300 mg) for the (HHL) reciprocal 
lattice plane and 19 single crystals (1 g) for the (HKO) 
zone. NiGa2S4 crystals were grown by the chemical vapor 
transport method followed by annealing in a sulphur at- 
mosphere. ICP analysis showed the stoichiometry to be 
within 1% of the nominal composition. The low-T lattice 
constants were a=6=3.624 A c=11.999 A a=/3=90° and 
7=120°. [3 

We first examine the development upon cooling of 
static SR spin correlations. Fig. [T] a) shows the T de- 
pendence of elastic scattering from powder and single 
crystals acquired with differing energy resolution. The 
single crystal and powder data on SPINS show a gradual 
onset of elastic scattering with decreasing T. For the pow- 
der sample, the increase in the intensity for 0.52<Q<1.0 
A"-'^ beyond that observed at T=40 K is a measure of 
the magnetic elastic scattering. The single crystal data 
were obtained through (HHl) scans fitted to a Lorentzian 
squared to obtain the integrated intensity (Fig. [2] 6)). 
On tightening the energy resolution from 5E=SQ ^eV 
(on SPINS) to 0.8 fieV (on HFBS), the characteristic T 
has decreased substantially and the onset of magnetic in- 
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FIG. 1: a) The T dependence of magnetic neutron scattering 
at Q=(0. 155, 0.155,1) for NiGa2S4 single crystals and powder 
taken on SPINS with E/=3.7 meV. h)-d) The magnetic spec- 
trum at 20 K, 10 K, and 1.4 K. The open and filled circles 
were taken with the crystals aligned in the (HHL) and (HKO) 
scattering planes and with the analyzer configured for 11° and 
5° acceptance respectively. The curves defined by the dashed 
line represents the energy resolution. 
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FIG. 2: a) - the elastic magnetic scattering versus momentum 
transfer along the L direction using 30 K as a background. 6) 
- the momentum dependence along the [110] direction with 
the vertical dashed line marking the commensurate (1/6,1/6) 
position, c) - the correlation length versus T measured in 2- 
axis mode. The horizontal bars are the resolution full-width. 



tensity upon cooling is sharper. The backscattering data 
indicate spin freezing on a timescale which exceeds r ^ 
fij^E ~ 1 ns. The strong correlation of the onset T with 
the energy resolution indicates that the elastic scatter- 
ing in NiGa2S4 is associated with rapid evolution with 
T of the time scale for incommensurate spin fluctuations 
rather than a conventional second order phase transition. 
These observations are reminiscent of observations at the 
spin glass transition in Cu-Mn. [14] The onset T for elas- 
tic magnetic scattering for 8E =0.8 /j,eV coincides with 
the broad maximum in magnetic specific heat (Cm) but 
is slightly above the 8.5 K susceptibility anomaly. Simi- 
lar observations were reported in strongly disordered 2D 
SCGO. For crystalline NiGa2S4 however, doping at 
the 1 % level was found to drastically suppress the low 
T maximum in Cm indicating an impurity concentration 
well below 1_% for our crystals, which do display the Cm 
anomaly. These results illustrate that our samples 
are very pure with minimal disorder present. 

We now discuss inelastic neutron scattering data that 
probe spin fluctuations versus T. The magnetic spec- 
trum at the critical wavevector Qo ~ (0.155,0.155,0) 
is displayed in Fig. [T] panels 6) through d) . The data 
were obtained on SPINS and corrected for a background 
measured at Q=(0.2,0.2,l) (open circles) and (0.25,0.5,0) 
(filled circles). At 20 K, there is no truely elastic scat- 
tering and the inelastic spectrum shows a broad peak 
in energy that can be described by a spin relaxational 
lineshape 5'(Q,Lj) cx [1/(1 - e- '^)]lu/{T^ + uj^) 



with ?ir=0.4(l) meV. At lower T, strong elastic (on the 
timescale of our measurement) magnetic scattering de- 
velops with two timescales; a resolution limited central 
peak and a broad relaxational response perhaps associ- 
ated with exchange or single ion anisotropy. The inelastic 
spectra indicate slowing spin fluctuations upon cooling, 
culminating in the development of a disordered staggered 
magnetization frozen on a time scale that exceeds ~ 1 ns. 

The momentum dependence of the elastic scattering 
is summarized in Fig. [2l Panels a) and b) illustrate 
correlations along the c-axis and within the a — b plane 
respectively. The solid line in panels a) and b) are fits 
to the following formula, which represents SR correla- 
tions within the a—b plane and correlations between NN 
triangular lattice planes. 

/(Q) = C^^^m2(l + 2acos(Q • c))e-<">'^' (1) 
2 2 (1 - sin^ cos^ (t)ql - cos^ 9qpA*('^/TT 

p{Q) is the Ni2+ form factor, £,= 26(3) A is the cor- 
relation length within the ab plane, (u^) is the mean 
squared displacement of Ni (taken to be 0.005 A^) [l7| , 
Qaf,=(H,K) is the in-plane wavevector, g is the Lande 
factor, (7^0)^ is 0.292 barn, A* is the area of the Bril- 
louin zone, and m = m(sin0cos(/)i: -I- e**^^ sin 9 sin (py + 
e^'I'^cosez) with X II to Qo=(0.155,0.155,0) and z \\ c. 
In terms of m, the time averaged spin on site r is given 
by (S(r)) = 2i?e{me'Q '"}. While S(r) depends on the 
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FIG. 3; Magnetic field dependence of the neutron scattering 
at Q=(0. 155,0. 155, L) measured on SPINS, a) - the field de- 
pendence of the change in peak intensity at L=0 and L=l. b) 
and c) - the intensity versus L at 10 and T with 25 K as a 
background. The horizontal bars are the resolution full-width. 



phase ^1 the magnetic scattering intensity for an unpo- 
larized beam (Eqn. [1]) does not. We assume the overall 
spin configuration (after domain averaging) is invariant 
under symmetry operations of the paramagnetic phase, 
eliminating effects of ip2- The parameter a=0.16(3) indi- 
cates nearest-neighbor inter-plane ferromagnetic correla- 
tions. The calibration constant C was determined from 
12 nuclear Bragg peaks ranging from (001) to (114) and 
indicates (5)^=0.26(5). The errorbar on (5)^ is due to 
uncertainty of the crystal mosaic. The strongly reduced 
sublattice magnetization (S—l) and is less than reported 
for a powder sample 11] and indicates strong quantum 
fiuctuations. 18| The dashed line in panel a) is a back- 
ground term cx resulting from non-cancellation of the 
Debye Waller factor affected elastic incoherent nuclear 
scattering. A Lorentzian squared was chosen to describe 
the momentum dependence of the in-plane correlations 
as it is normalizable in 2D and has been used to describe 
SR correlations in random field magnets. 19] 

Eqn. [T] convolved with the resolution fits the complete 
L dependence. The broad maxima for integer L indicate 
ferromagnetic correlations between NN planes. Two spin 
configurations fit the single crystal data equally well. An 
out of plane spiral with </> = 0° and 6 — 57(5)° and an in- 
plane spiral with (p — 41(5)° and 9 = 90°. Higher angle 
peaks in the powder data are better described with spins 
in the basal plane. jll| 

We now examine in-plane spin correlations. Fig. [2] b) 
shows a scan along the (HHl) direction. Instead of the 
conventional 120° structure, the data are more closely 
described by a 60° structure as expected for dominant 
third NN interactions (J3) consistent with first princi- 
ples calculations. [20| In addition, the peak is clearly 
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FIG. 4: The momentum dependence of magnetic neutron 
scattering measured on SPINS, a) is the dispersion and b) 
-d) are constant energy scans. Data at ?j,tj=0.5, 1.0, 2.0, and 
3.0 meV were taken with E/=3.7 meV and hLj—2.5 and 3.5 
meV with E/=5.1 meV. The analyzer was configured for 5° 
horizontal acceptance. The horizontal bars are the resolution 
full-width. 



shifted below the Q=(l/6,l/6) commensurate position. 
As detailed below, the incommensurability may result 
from competing first and third NN interactions. 

The T dependence of the inverse in-plane ^ is plot- 
ted in Fig. [2]c) as obtained on the BT2 spectrometer in 
two-axis mode with Ei=14.7 mcV. The data represent an 
approximate measurement of the instantaneous correla- 
tion length. ^ increases upon cooling saturating at 26(3) 
A with no finite-T anomaly as would be expected at a 
conventional second order phase transition. 

To better understand the origin of the SR order, we ex- 
amined the spin response to an in-plane magnetic field. 
Fig. [3] b) and c) show scans along c* at 10 T and T. 
With increasing field, there is a reduction in the ampli- 
tude of the intensity modulation indicating a decrease in 
the c axis correlations (reflected by a in Eq. 1). Fig. [3]a) 
shows the field dependence of the intensity at two points 
along the L-scan: L=0 where the intensity increases and 
L=:l where it decreases with field. Quantitatively, the fit 
to Eq. [T] represented by the solid lines in panels b) and c) 
give a(0r)/Q;(10T)=1.5(2). Therefore, inter-plane cor- 
relations are suppressed in a magnetic field. In substi- 
tutionally disordered magnets, the application of a mag- 
netic field produces a random field which couples linearly 
to the order parameter. 
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the high field data is that in-plane pinning is strength- 
ened in a field at the expense of inter-plane correlations. 

We now examine magnetic excitations in the low-T 
limit from the SR ordered frozen state at T=1.5 K. Figs, 
m (b)-4(d) display constant energy cuts along the (HHl) 
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direction for energies 0.5 meV, 2 meV and 3.5 meV. The 
inelastic spectrum is summarized in Fig. |4]in the (HKO) 
scattering plane. For increasing energy transfer the cor- 
related magnetic scattering broadens in momentum and 
then disappears at higher energy transfers. The data can 
be described by two peaks displaced symmetrically from 
Q=(0. 155,0. 155) and separating in proportion to energy 
corresponding to a velocity ?ic=29(3) meV A. Fig. |4]a) 
shows the peak positions at each energy. Owing to the 
short-range correlations observed, the peaks are signifi- 
cantly broader than resolution. 

The spin-wave velocity, Qcw and S can be described 
in terms of a weak NN ferromagnetic coupling (Ji) and 
a dominant third NN coupling (J3). The peak position 
Qo is associated with the minimum in the exchange field 
defined as Jq = "^^^ J„e'^~*'i ''"^ where J„ is the exchange 
constant for a neighoring spin located at a distance r. 
The spin- wave dispersion can be obtained from [22 1 



change. This model, however, contrasts with the lack 
of long-range spin order. Quantum fluctuations, and/or 
anomalous domain dynamics 0, with dilute strong 
pinning [lo| are avenues that should be further explored 
to determine why long range order does not occur in high 
quality single crystalline NiGa2S4. Irrespective of the 
eventual resolution of that puzzle, our results indicate 
that the S=l antferromagnetic on a triangular lattice is 
a marginal case and the collective properties cannot be 
understood in terms of conventional classical theory as 
in large- S systems. 
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q is defined as the displacement in momentum away from 
the ordering wavevector at Qo and (S) sublattice magne- 
tization. To reproduce both Qo=(0. 155(3), 0.155(3)), the 
spin- wave velocity ?ic=29(3) meV A, and ecw=80{2) 
we find J3=2.8(6) meV and Ji=-0.35(9) J3. These val- 
ues result in the curve in Fig. IHa), give Qcw — -87 K 
and Qo=(0. 153, 0.153), consistent with experiment. The 
weakly ferromagnetic Ji is compatible with the cuprates 
for Cu-O-Cu bond angles less than - 100°. [H, Hj] The 
NN Ni-S-Ni bond angle in NiGa2S4 is - 97°. 

Quasi-2D gapless spin waves are also inferred from heat 
capacity measurements. [ll| In the low T limit C/R — 
(i/n33/2a^fc|/(27rc2?i^))f (3)r2, where a is the lattice con- 
stant, ks is the Boltzmann constant, and ^(3) ^ 1.202. 
From the low-T Cm we obtain he =5.3 ^/nU meV A. 
Assuming, in any given volume element, there are only 2 
magnetic critical wavevectors in the Brillouin zone {n=2) 
and the peak in Fig. 1 (d) reflects easy plane anisotropy 
(^=1) we find he ^ 8 meV A. The discrepancy between 
Cm and neutron data could result from a range of veloci- 
ties in the disordered state or a slow mode that does not 



contribute to neutron scattering. [25|, l26| 

We have shown that NiGa2S4 is a 2D triangular lat- 
tice antiferromagnetic which does not develop long range 
order down to T/Qcw ~ 2% even in high quality sin- 
gle crystals. Instead spin-freezing occurs for T~ 10 K 
and the low T state consists of 2D SR incommensurate 
order. Many of the physical properties including the in- 
commensurability, the spin-wave velocity and the Curie- 
Weiss temperature can be accounted for by strong third 
NN antiferromagnetic and weak NN ferromagnetic ex- 
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